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Abstract 

Motivated by rolling adhesion of white blood cells in the vasculature, we study how cells move 
in linear shear flow above a wall to which they can adhere via specific receptor-ligand bonds. Our 
computer simulations are based on a Langevin equation accounting for hydrodynamic interactions, 
thermal fluctuations and adhesive interactions. In contrast to earlier approaches, our model not 
only includes stochastic rules for the formation and rupture of bonds, but also fully resolves both 
receptor and ligand positions. We identify five different dynamic states of motion in regard to 
the translational and angular velocities of the cell. The transitions between the different states 
are mapped out in a dynamic state diagram as a function of the rates for bond formation and 
rupture. For example, as the cell starts to adhere under the action of bonds, its translational and 
angular velocities become synchronized and the dynamic state changes from slipping to rolling. 
We also investigate the effect of non-molecular parameters. In particular, we find that an increase 
in viscosity of the medium leads to a characteristic expansion of the region of stable rolling to the 
expense of the region of firm adhesion, but not to the expense of the regions of free or transient 
motion. Our results can be used in an inverse approach to determine single bond parameters from 
flow chamber data on rolling adhesion. 

PACS numbers: 82.39.-k, 83.10.Mj, 87.17.Aa 
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I. INTRODUCTION 



Rolling adhesion of white blood cells (leukocytes) along the vascular endothelium plays 
a key role in the immune response and is a prominent example for the interplay between 
transport and specific adhesion in biological systems . Random encounters between 

cell and vessel wall lead to the formation of initial bonds based on adhesion receptors from 
the selectin-family. Because selectins are characterized by fast association and dissociation 
kinetics, new bonds are readily formed on the downstream side and old bonds are contin- 
uously ruptured at the upstream side. If the processes of bond formation and rupture are 
sufficiently balanced, rolling adhesion results. As rolling velocity is smaller than the velocity 
of a cell moving freely in hydrodynamic flow, this mechanism allows the leukocytes to more 
efficiently survey the vessel wall for appropriate molecular signals. The main signals in this 
context are chemokine molecules, which indicate the presence of an infection and lead to 
firm adhesion based on the adhesion receptors from the integrin-family. Firm adhesion in 
turn is a prerequisite for transendothelial migration into the surrounding tissue. Similar 
mechanisms are used by stem and cancer cells to disseminate in the body through the blood 
flow. 

The main tool for investigating rolling adhesion under controlled conditions are flow 
chambers |4i]. There receptor-carrying cells suspended in an aqueous solution are perfused 
in linear shear flow above a ligand-coated wall. For example, it was shown with flow cham- 
ber experiments that selectin bonds have fast kinetic rates and that the dissociation rate 
increases with force (slip bond), thus explaining their superior function in mediating rolling 
adhesion . Recently it was demonstrated by a combination of atomic force microscopy 
and flow chamber experiments that the lifetime of single selectin bonds shows a biphasic 



response under force 
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I . Such a catch-slip bond behavior might have evolved to avoid ad- 
hesion of leukocytes under static conditions. Indeed it has been found early in flow chamber 
experiments that leukocytes adhere only above a critical threshold of shear [ol. In addi- 
tion to speciflc properties of the molecular bonds, transport processes might also play an 



important role in creating the shear threshold 10|, lUj. Both rolling adhesion 12| and the 



shear threshold 



131 ] have also been demonstrated in cell-free flow chamber experiments with 



ligand-coated micron-sized beads. Although cell-free rolling is more erratic than leukocyte 
rolling, this indicates that cellular features (like cell deformability) are not essential to rolling 
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adhesion. Flow chamber experiments also allow to study the exact effect of non-molecular 
parameters, for example of changes in the viscosity of the medium. By adding an inert sub- 
stance like the sugar ficoll, it is possible to change shear stress but not shear rate in a flow 
chamber experiment, thus dissecting the respective roles of force and transport 1^, [ll|, 1^ . 

The standard observable in flow chamber experiments is the translational velocity as a 
function of time 0, [is], 1^. From this time series further variables can be derived, for 
example average and standard deviation of velocity. In the case of low ligand density, a 
stopped cell is most likely held by a single bond and detaches without rebinding. Then the 
time series can be used to measure single bond lifetime. At higher ligand density, rolling 
occurs. Because rolling is never smooth but erratic due to the stochastic processes on the 
molecular level, then the time series can be used to measure pause time distributions. The 
velocity time series can also be used to define different states of motion of the cell. Typically 
a cell is said to undergo rolling whenever its mean velocity (averaged over some seconds) 
significantly decreases compared to the free (hydrodynamic) velocity. If no motion can be 
detected for several seconds, the cell is considered to be in firm adhesion [l^. 

Because rolling adhesion is not only of large physiological importance, but also character- 
ized by an intricate interplay of different physical factors, it has long been subject to intense 
modeling efforts. If molecular effects are of interest, then such modeling efforts typically 
start with physical models for bond association and dissociation [l^]. Combined with the 
hydrodynamics of a sphere in front of a wall, they lead to algorithms known as adhesive dy- 
namics 13]. This approach has been applied before to different aspects of rolling adhesion, 
e.g., the interplay of two receptor systems 19 1 or the effect of catch bonds |2o|. Recently 
we have developed a new variant of this algorithm which in contrast to earlier approaches 
fully resolves the spatial positions of the receptors on the sphere and the ligands on the wall. 
Using this approach, we were able to predict the efficiency of initiating cell adhesion in shear 



flow as a function of the density and geometry of the receptor and ligand patches 



large modeling effort has also been spent on the role of cell 
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deformability 23 



24 



21 
2 



22|. A 

m 



27| 



291]. Here we focus on the case of 



and the interaction between multiple particles |27l . 
moderate shear flow and small numbers of cells, when cell deformability and hydrodynamic 
interactions between cells are not relevant {2^ . 

A convenient way to present the results from adhesive dynamics simulations of rolling 
adhesion is the calculation of state diagrams which predict different types of motion over 
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a large range of model parameters 30|, ImI- Because in flow chamber experiments one 



usually measures only translational velocity as a function of time, these state diagrams have 
been determined before based only on the translational velocity simulated as a function 
of different molecular parameters. However, computer simulations also allow to track the 
angular velocity of a cell, thus doubling the number of degrees of freedom that could be 
compared. For example, in physical terms cell rolling means that translation and rotation 
are synchronized. Although this fact was already noted in the pioneering paper on adhesive 



dynamics [18j, it has not been systematically exploited due to the lack of experimental 
data. Here we present a detailed analysis of the different dynamic states of rolling adhesion 
which is based on the simultaneous tracking of both translational and rotational degrees of 
freedom. In contrast to earlier work in this fleld, we fully resolve the spatial positions of 
receptors and ligands and calculate the state diagrams as a function of the on- and off-rate. 
Moreover, we explore the role of different external parameters, including the viscosity of the 
medium. 

The paper is organized as follows. In Sec. [Ill we explain our numerical method, which 
is a combination of a Langevin equation accounting for hydro dynamic and thermal forces 



in the limit of small Reynolds numbers (Stokesian dynamics) 



32| and adhesive dynamics 



accoun ting for the forces caused by the formation of adhesive bonds between cell and sub- 



strate 



18| . By taking care of the diffusive motion of the sphere, we are able to explicitly 



resolve both receptor and ligand positions. As a first application of the hydrodynamic part 
of our model, we explain the physical difference between slipping and rolling. We close this 
section with a detailed description of the parameters relevant for the algorithm. In Sec. IIIII 
we point out, using a simple analytical description, how the cell slows down and simultane- 
ously synchronizes its translational and angular velocity under the action of bonds. Then 
we discuss the mean velocities of the cell as a function of the on- and off-rate of the receptor- 
ligand bonds. Based on this we define five distinct states of stationary motion for a cell in 
a fiow chamber. In Sec. [IV] we identify these different states in state diagrams displaying 
the dependence on the internal bond parameters, i.e., the on- and off-rate. In addition, the 
effect of external parameters to the location of the states in the space of different on- and 
off-rates is discussed. In the closing section. Sec. |Vl we argue that molecular parameters 
can be much better extracted from flow chamber experiments if both translational and ro- 
tational degrees of freedom are measured, as it is done in our simulations. Furthermore we 
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FIG. 1: Model system and adhesive dynamics. A rigid sphere of radius R moves above a wah. 
The shear flow is linear with shear rate 7. Receptors are modeled as sticky spheres with capture 
radius vq <^ R on the surface of the sphere. Ligands are modeled as dots on the boundary wall. 
A receptor-ligand bond forms with rate kon whenever a receptor patch has some overlap with a 
ligand. Bonds are modeled as Hookean springs and exert a force on the sphere that is proportional 
to the distance between the receptor position and the ligand position r/. They rupture with a 
force-dependent rate koff{F). 

propose experiments based on recent nanotechnological developments that could result in 
flow chamber data for the angular velocities of cells or microspheres. 

II. MODEL AND SIMULATION METHOD 
A. Stokesian dynamics 

As a simple model system for a cell or a microsphere in a flow chamber we use a rigid 
sphere of radius R that moves above a wall (cf. Fig. [1]). For objects as small as white blood 
cells (with diameters < 10 fim) the typical value for the Reynolds number is much less than 
one and inertia can be neglected (overdamped motion). Therefore, the flow around the cell 
is laminar and well described by the linear Stokes equation. Cells are usually observed only 
close to the wall, that is at a distance to the wall which is much smaller than the separation 
between the two parallel glass plates used in a flow chamber. In this region the unperturbed 
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flow profile is approximately linear. Throughout this paper we therefore consider the sphere 
to flow in simple shear flow, i. e., we consider the shear rate 7 to be a constant. 

In the situation under consideration the sphere is not only driven by the imposed shear 
flow but also by direct forces exerted on the cell. These include a constant gravitational drift 
force which results from a slight density difference between the cells and the surrounding 
medium. For flow chamber experiments, this implies that most cells adhere to the bottom 
plate, which usually is also the plate coated with ligands. In our context, the strongest 
forces result after initial binding from receptor-ligand bonds that pull on the sphere. A 
third source of force are thermal forces which are ubiquitous for objects in the micrometer 
range and lead to diffusive motion. This kind of motion can be described by an appropriate 
Langevin equation. It is convenient to introduce a six-dimensional state vector X which 
includes both the three translational and the three rotational degrees of freedom. The first 
three components of X denote the Cartesian coordinates with respect to a flow chamber 
based coordinate system (cf. Fig. [1]). The last three components describe the rotation that 
maps a coordinate system fixed to the center and the orientation of the sphere to the flow 
chamber fixed coordinate system. Similar compact notations are introduced for velocities 
and forces. The symbol U denotes a six-dimensional velocity with the first three components 
being translational velocities and the last three being angular velocities. F is the combined 
force and torque vector. Using this notation, the Langevin equation describing the motion 



of the cell reads 22 1 



X = U°° + M(F^ + F^) + A;BTVM + g/. (1) 

Here, M is the 6x6 mobility matrix. U°° oc 7 is the velocity of the unperturbed shear flow 
at the position of the center of the sphere. F*^ oc 7 is the shear force, which results from 
the hydrodynamic interaction between the cell and the wall. F^ denotes all direct forces 
(torques) acting on the sphere, like gravity and bond forces. The terms /c^TVIVI and 
describe the effect of thermal noise, gf is Gaussian white noise with 

id) = 0, {gkl) = 2kBTaM6{t - 1') . (2) 

Here Ta is the ambient temperature and ks Boltzmann's constant. The mobility matrix for 
a sphere above a wall depends the distance between sphere and wall [sS^ . Thus, the noise is 
multiplicative, leading to the gradient term VM in Eq. ([1]). We interpret the noise gj in the 



usual Stratonovich sense. However, Eq. ([T]) is written in its Ito version [3J, l35j, as indicated 
by the superscript /. This allows to directly derive a simple Euler algorithm for the update 
step AX of the configuration of the sphere during a time step At. In non-dimensional form, 
using the radius R of the sphere as the length scale, the inverse shear rate I/7 as the time 
scale, and GnriR'^'y as the force scale, the first order discretized version of Eq. ([T]) reads 

AXt = (U°° + M(F^ + F))At + -^VMAt + ^J^g{At) + OiAf). (3) 

Here the dimensionless Peclet number Pe = {GrcqR^'^) / (ksTa) describes the relative impor- 
tance of convective versus diffusive motion of the sphere. The first two moments for the 
Gaussian white noise now read 

(g(At)) = 0, (g(At)g(At)) = 2MAt . (4) 

The algorithm Eq. (j3l) is also known as Stokesian dynamics and has been derived, e.g., by 



Brady and Bossis [321] and for vanishing shear flow by Ermak and McCammon [36|]. A more 
detailed description of our algorithm is given in Refs. [s^] . In order to obtain accurate 
results when simulating the motion of the sphere it is essential to properly calculate the 



mobility matrix M and the shear force F'^. In our simulations we use a numerical method 



described by Jones et al. [33|, |38j which allows to accurately calculate the components of 
the mobility matrix for a single sphere above a wall for arbitrary sphere wall distances. The 
values obtained by this method agree very well with the classical results given for some 



tabulated height values by Goldman et al. [39 1. 



From Eq. ([3]) one readily sees that diffusive motion is less relevant if the Peclet number 
Pe is large. Using typical values for white blood cells (cf. Tab. [T]) one finds that Pe is 
of the order 10^. Therefore, regarding the motion in flow direction thermal motion can be 
neglected on the scale of the cell radius. However, for the motion in z-direction the only 
force on a freely flowing cell is gravitation, which is small due the small density difference 
between the cell and the surrounding medium and therefore comparable to thermal forces 
in this direction 22|. For this reason, in our work we use the full Langevin equation from 
Eq. ([1]). Another reason why we keep the diffusion terms in Eq. ([3]) despite the large Pe 
values has to do with the size of the receptors on the surface of the cell and will become 
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Z-R 

FIG. 2: Velocities of the sphere in the hmit of deterministic motion, Pe oo. Translational 
and angular velocity U and respectively, and their ratio RQ/U are shown as a function of the 
separation z — R between sphere and wall. The thin lines depict the asymptotic behavior for z ^ R 
and z ^ oo, respectively. U and are plotted in units of R'j and 7, respectively, z is measured in 
units of R. 

more apparent in the next subsection, Sec. Ill B[ 

For the following it is important to consider how the translational velocity U (in x- 
direction) and the angular velocity fl (for rotations about the y-axis) behave in the limits 
of small and large separation from the wall. If the translational and angular velocities were 
completely synchronized, the ratio RQ/U should be equal to unity. This would correspond 
to rolling in a macroscopic sense, e.g. for a sticky sphere rolling down an inclined plane or 
for a car wheel on the street. We now show that the situation is different for a cell in free 
hydrodynamic motion. For z denoting the height of the center of the sphere above the wall, 
z — R denotes the gap between sphere and wall. In Fig. [2] we plot U and Q a.s a. function 
of 2; — i? in the limit of deterministic motion, Pe 00. For the sphere far away from the 
wall, z — R 00, we have the case of free flow, thus U = 'yz and Q = 1/2. Therefore the 
ratio RQ/U vanishes in this limit. For the sphere approaching the wall, z — R —>■ 0, both U 
and Q slowly approach zero due to the hydrodynamic no-slip boundary condition. The ratio 
Rfl/U however approaches a finite limit RQ/U ~ 0.5676 as computed by Goldman et al. 
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39l |. Together these results show that hydro dynamic interactions increasingly synchronize 
translational and angular velocities as the cell approaches the wall. However, the ratio 
Rfl/U never reaches the value 1 as it would for rolling in a macroscopic sense. Therefore, 
the hydrodynamic coupling between the cell and the wall is not strong enough to lead to 
rolling and a freely moving cell is always slipping. 

B. Bond dynamics 

We now include receptors, ligands and receptor-ligand bonds into our model. Throughout 
this paper we call the cellular part of a bond receptor and its counter part on the wall ligand. 
Before a cell receptor can react with a wall ligand to form a complex (bond), there must be 
a physical transport process which brings the two components to close proximity. Formally 
bond formation can be separated into a transport and a reaction step using the notion of an 
encounter complex . The encounter complex is formed whenever a receptor and a ligand 
are less than the capture radius tq away from each other. Therefore, we model receptors as 
reactive patches on the sphere surface having a spherical capture range of radius tq. The 
capture length tq bridges the gap between the continuum approach followed by implementing 
linear hydrodynamics and the discrete nature of the receptor molecules. Ligands are modeled 
as dots on the boundary wall. Because we explicitly resolve receptor and ligand positions, 
it is essential to also include diffusive motion in the algorithm, Eq. (JSj). Although diffusion 
plays a minor role for the relative position of the cell, it affects the positions of its surface 
receptors, which are of much smaller size than the sphere itself. 

During the time an encounter complex exists it can react to the final receptor-ligand bond 
complex with the on-rate kon- Reversely, any bond complex can rupture into the encounter 
state with the off-rate kofj. Experimentally it has been found that the dissociation rate 
depends on the physical force acting on the bond complex P, 0]. These early experiments 
agreed nicely with the simplest model for bond dissociation under force, which had been 



conceived first by Bell 17|: 



koff{F) = k^eMF/Fd), (5) 
with fco the dissociation rate at zero pulhng force, F the force on the bond and Fd the 
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detachment force scale. The Bell model, Eq. ([5]), can be rationalized using Kramers theory 
as being a thermally activated escape over a transition state barrier in the presence of an 
external force 41, 42, 43, 4^. Although recent evidence suggests that the dissociation 
rate shows a more complicated force dependence at small forces 3, Eq- O has been 
demonstrated to properly describe the dissociation process for selectin bonds in the high 
force regime [ij]. Here we use the Bell model for conceptual and computational simplicity. 

In order to include the probabilistic nature of bond formation and rupture, the algorithm 
Eq. ([3]) has to be extended to include rules that take care of these processes. Such rules 
have been first set-up by Hammer and co-workers [18| and been refined many times to 
model various aspects of leukocyte rolling. They are now known as adhesive dynamics. 
Here, we briefly explain the main idea behind these rules. A detailed description of our 
implementation of these rules is then given in the appendix. Sec. \M In order to simulate 
the motion of a cell under the action of bonds, the configuration of the sphere is updated 
at each time step At according to Eq. ([3]). If, after some update, a receptor patch on the 
surface of the sphere overlaps with a wall ligand a receptor-ligand bond forms with rate kon, 
i.e., the probability for bond formation during time step At is 1 — exp{—konAt). Existing 
bonds are modeled as harmonic springs, i.e., the force along these bonds is proportional 
to the bond extension, essentially given by ||r,. — r;||, with r,. and r; being the receptor 
and ligand position, respectively (see Fig. [1]). Force and torque resulting from extended 
bonds enter then the configuration update equation, Eq. ([3]), via the direct force term F^. 
Similarly to bond formation, the probability for bond rupture during time step At is given by 
1 — exp{—koffAt), where fco// is the force dependent dissociation rate according to Eq. (JSj). 

In contrast to earlier implementations of adhesive dynamics, we explicitly resolve both 
receptors and ligands in space, rather than considering a wall homogeneously coated with 
ligands at constant density. One immediate advantage of our method is that it avoids a flow 
rate-dependent rate of bond formation 451]. In order to be able to spatially resolve both 
ligands and receptors, it is necessary to include the Brownian motion of the cell. If one 
considered only the deterministic part of Eq. ([T]), at low densities of receptors and ligands 
it could happen that a receptor never finds a ligand. Although these extensions of adhesive 
dynamics lead to increased computational effort, they are closer to real biological systems, 
which on the molecular level are in permanent thermal motion. In particular, in the future 
our implementation will allow us to model receptor-ligand kinetics in more molecular detail. 
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As will be discussed later, our algorithm also opens up the perspective to compare adhesive 
dynamics simulations to flow chamber experiments using substrates with nanopatterned 
ligand. 



C. Parameters 

Our model contains thirteen different dimensional parameters. With i?, I/7, GiTTjR'^'y 
being the natural scales of length, time and force, respectively, we are left with ten nu- 
merical (dimensionless) parameters appearing in the algorithm. Typical values for these 
(both the dimensional and dimensionless) parameters are listed in Tab. [H The parameters 
i?, Ta, 7, r], Ap influence the flow properties and besides viscosity we keep these parameters 
flxed. For the ambient temperature we choose room temperature = 293 K. For flow 
chamber experiments with cells, the choice Ta = 310 K would be more appropriate, but the 
exact choice of this value, which affects bond dissociation kinetics and effective strength of 
diffusion, turns out to be irrelevant for the physiologically relevant parameter values chosen 
here. For the Stokes radius R of the cells we use i? = 4.5 /im which is about the measured 



radius of neutrophils, a main type of leukocytes undergoing rolling adhesion [46|. The Nr 
receptor patches are randomly distributed on the cell surface. The receptor patches might 
be identifled with cell microvilli, which are membrane projections to whose tips the selectin 



receptors are localized 5J]. The number of microvilli on a leukocyte varies from several 



hundreds [52|] up to 10,000 [18|, l5l|- In our simulations, we use A^,. in the order of 10^. Al- 
though usually several receptors can be found on the microvilli tips, we allow only one bond 
per receptor patch, i.e., Nr is also the total number of receptors. Receptors and ligands 
are spatially extended in the nm-range and their binding sites diffuse within some region 
around their linkage. The diffusion constant of a nm-sized object is about three orders of 
magnitude larger than that of the cell itself. Therefore, it is sufficient to account for the 
spatial distribution of the location of binding by introducing a capture sphere with capture 



radius tq = 50 nm, which is about the combined length of ligand and receptor [53|, |55|] . In 
regard to the ligands, we consider them being distributed on a square lattice with lattice 
constant d. Here, d is obtained from d = a/ 1 / A"; , with A"; being the average number of 
igands per /im^. In flow chamber experiments A^; typically varies between (1 — 100)//im^ 

3- 
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Parameter — typical value meaning reference 
non-dimensionalized (dimensionless) 

R ^ 1 4.5 ... 5 /im radius [6, 46] 

7^ 1 50 ... 150 Hz shear rate [14] 

T(j 293 . . . 310 K ambient temperature 

77 1 ... 3 Pa s viscosity [14] 

Ap 50 kg/m^ density difference [47] 

«; 10"^ . . . 10"^ N/m bond spring constant [18, 48, 49] 

-^k/6717]R^ {m-\..io^) 

kon 10^.. 10^ Hz on-rate [10] 

vr = fc„„/7 (10-3. . . 10) 

ko 0.5 . . . 300 Hz unstressed off-rate [6, 50\ 

^eo = V7 (10-4... 103) 

To 50 nm capture radius 

^ ro/R (10-2) 

c? 0.1 ... 1 /im ligand-ligand distance [5, 6] 

^d/i? (0.02... 0.2) 

Xc 2 ... 4 ■ 10-^^ m reactive compliance [6] 

-^GnrjR'^'yXc/kBTa (0.1... 0.6) 

Nr 50 . . . 5000 Number of receptors [18, 51, 52] 

/imin 15 nm minimum cell height [53] 

^ /^m»n/i^ (3-10-3) 



TABLE I: Parameters used for the adhesive dynamics simulations. If no extra symbol for the di- 
mensionless quantity is defined we use the same symbol for both the dimensional and dimensionless 
representation of this quantity. 

The on-rate kon for single bond formation from a receptor-ligand encounter in Tab. [T] 
is given in units of Hz. Experimentally, it is very difficult to determine this rate directly. 



Using three-dimensional affinity data, an upper limit has been estimated to be 10^ Hz 10|. 
For the force dependent off-rate fco// we have the two Bell parameters k^^Fd, where the 
detachment force is Fd = ksTa/xc with reactive compliance Xc- Both, unstressed off-rate 
fco and reactive compliance Xc have been measured for different selectin bonds The 
reactive compliance for L-selectin bonds is = 2 - 10-^^ m [6]. This corresponds to a typical 
detachment force of 200 pN. 

For closed bonds we use the linear force extension curve explained in Sec. Ill BI with spring 



constant n. For the P-selectin-PSGL-1 complex Fritz et al. 48] measured a value of hrl 



5.3 ■ 10 3 N/m. Recently, also the role of microvilli elasticity was discussed 3l|. Shao et al. 
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(a) (b) 

FIG. 3: Snapshots of leukocyte motion, (a) Velocity of a leukocyte (neutrophil) rolling on L-selectin 
ligand PNAd with 60 sites per fim^. (Data obtained from Ref. Fig. lA]). (b) Translational and 
angular velocity (U and RQ, respectively) of a sphere with radius R = 5 /xm measured in adhesive 
dynamics simulation at the same ligand density (with on-rate kon = 60 Hz and unstressed off-rate 
ko = 6.8 Hz). 



determined the spring constant of microvilli in the low force regime to be = 4.3-10-^ N/m 
491]. The total spring constant k = k„ivI^rlI {i^mv + i^rl) of the series of microvilli and bond 
would then be dominated by the microvilli spring constant. Besides the bond forces, we 
include only the buoyant force due to the small density difference Ap between the cell and 
the surrounding medium. Other nonspecific repulsion forces arising from electrostatic and 
steric stabilization forces are effectively taken into account by introducing a simulation rule 
that the cell can approach the wall only up to a distance of /imm = 15 nm 53|]. In a 
physiological system, this would correspond to the thickness of the glycocalyx, a protective 
layer of sugar covering the endothelium in blood vessels. 



III. DEFINITION OF THE DYNAMIC STATES OF MOTION 



A. Experiment and simulation 

Having completed the model definition, we now start to analyze the simulation results. 
For the following we first have to clarify what is meant by the velocity of the cell. As we 
include Brownian motion, the velocity of the cell U{t) cannot be its instantaneous velocity, 
because for the trajectory of a Brownian particle \\m5t^Q{X{t + 6t) — X{t))/5t is not a well- 



defined quantity |56|. Thus, we define the velocity of the sphere at time t by a difference 
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quotient U{t) := {X {t + At) — X {t)) / At with some time interval At. Throughout this paper, 
we choose At = 0.02 s, which corresponds to a camera resolution of 50 Hz. The angular 
velocity for rotations about the y-axis Q{t) is defined in a similar way. 

In Fig. [3] we compare in a qualitative way the data obtained in a numerical simulation 
of leukocyte motion with data of a rolling cell obtained in a flow chamber experiment. In 
Fig. [3^ the translational velocity U (denoting the velocity in the direction of imposed shear 
flow) of a rolling leukocyte for some period of time is shown as experimentally measured by 
Alon and co-workers The rolling state is identified by a strong decrease of the average 
cell velocity compared to the average velocity of a cell moving freely in hydrodynamic flow. 
Fig. [3)3 shows a representative trajectory from our simulations. No attempt has been made 
to fit this data set to the experimental one. Yet it is clear that both data sets show very 
similar features, including the strong decrease of velocity upon binding and the subsequent 
bursts in velocity, which correspond to stochastic rupture events at the trailing edge allowing 
for cell movement. In contrast to the experimental data, the simulation data also records 
the angular velocity Rfl. The simulation reveals that the RQ and U curves collapse onto 
one curve as soon as the cell binds for the first time (at t ^ 1.2 s). Before initial binding, the 
cell slips over the substrate with RVt/U < 0.57 as discussed in Sec. Ill A[ This observation 
motivates us to define the dynamic state of rolling by RQ/U ^ 1, in accordance with the 
common understanding of this term in macroscopic mechanics. However, before we develop 
this idea in more detail in Sec. IIII CI we first show with an analytical calculation how the 
action of a bond synchronizes translational and rotational velocities. 

B. The stopping-process: emergence of rolling 

In order to understand how the cell comes to a stop after the first bond has been formed, 
in the following we consider a set of simplified equations of motion. We reduce our analysis 
to a two dimensional case in which the cell is only allowed to translate in the x-direction and 
to rotate about the y-axis. These two degrees of freedom are called x for the translational 
and 6 for the rotational degree of freedom, respectively. In regard to the z-direction we 
assume that the cell moves at constant height above the wall {z = 0) with z = R + rQ. 
Furthermore, we neglect Brownian motion, i.e., we consider the limit Pe —* oo. Then, the 
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FIG. 4: A cell is stopped by a tether force F^. r/ and r,. are ligand and receptors position, 
respectively. x,6 define two degrees of freedom, x the angle between the bond and the wall. In 
equilibrium (i. e. U = x = 0,O, = 9 = 0) the shear force and the shear torque are balanced 
by the tether force and its respective torque. Not shown are gravitational force and repulsive forces 
from the substrate that compensate all downward acting forces. 



equations of motion for the two coordinates read (with x = U and 9 = Q, cf. Sec. Ill All 



X \ I Fb,x \ / Uhd 

d I \ TB,y I \ ^hd 



(6) 



with Fb,x the x-component of the bond force and Tb ^ the y-component of the torque which 
is due to the bond (see Fig. H]). Because the bond is modeled as a linear spring, it has 
mechanical properties and we call it a tether. The matrix M^g denotes the (x^)-sector of 
the mobility matrix introduced in Eq. ([1]). If no bond is formed, then Fb^x = Ts^y = 0, and 
U = Uhd and Q = Qhd- The free velocity Um is often referred to as the hydrodynamic velocity 
of the cell. From Fig. H] we read off tether force and tether torque and T^, respectively: 



FBit) = K 



^ -x{t) + Rsm{e{t))^ 


^ Rcos{9{t))-z J 



^ -Rsm{e{t))^ 



TBit) = rxFB, f 



V" 





-i?COS 



(7) 



In the following we distinguish two different cases with respect to the resting length of the 
bonds. We start with the case were we assume the resting length to be zero. With the cell 
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moving at a constant height z = R + tq, at t = aheady a small bond force in 2;-direction 
exists (but no torque). Briefly after bond formation the cell has not moved significantly and 
x/R,9<ti 1 holds. Thus, the quantities FB,x,TB,y can be approximated as 

Fb,x = f^{RO -x) + 0{e^), TB,y = k{Rx - zRd) + 0{6^) + O^xO'^). (8) 

Reinserting these approximate expressions into Eq. ([6]), we obtain a first order linear differ- 
ential equation for x, 



C := M,,X, X := « , (9) 






which can readily be solved with the proper boundary conditions (a;(0) = ^(0) = 0). The 
determinant of the matrix C is det C = k,'^R{z — R) det M^q. From the dissipative nature of 
the mobility matrix expressed, e.g., in Eq. ([2]), it follows that M^e is positive definite, which 
is then also true for the matrix C. As the diagonal elements of M^g are significantly larger 



A+-(c)ii / A_ - A+ ' ' A_-(c)n / A_ - A_ 

(C)l2 / \ (C)l 



than the off-diagonal element [33|, |38|], one easily finds that trC is negative. Therefore, C has 
two negative eigenvalues, which we denote in the following by A±. As z — R R one can 
approximate A+ ~ trC and A_ ~ det C/trC. These two eigenvalues represent two different 
timescales with |A_| ^ |A+|. With this the solution to Eq. ([9]) can be written as 

where the asymptotic solution after the cell has stopped is given by 



(11) 



Eq. (fTTj) is the linearized version of the force and torque balance condition at mechanical 
equilibrium. In non-linearized form the force and torque balance equation reads 



1 -^oo 


1 - -(^'-' 1 
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Then, for z = IMR one gets Fb,x ~ -1-7 ■ GirviR^ and TB,y ^ -0.6 ■ 67rr]^R^ |39|. The 
relation to the quantities Xoo, X is given in Ref. [gI 

||Fb||cosx= ilFsll, ||FB||xooSinx= IIT5II +i?||Fs||. (12) 

The relation with ^^oo is given by the purely geometrical relation Xoo/R = sin^oo + (1 — 
cos ^00)/ tanx (cf. Fig. Hj). In contrast to the linear version, Eq. (fT2|) cannot be solved for 

The initial velocities of the cell are the free hydrodynamic velocities with R6{0) /i(0) ~ 0.5 
so that X and 6 do not increase at the same speed (see Fig.Ek). Shortly after bond formation, 
t <^ 1/|A+|, exp(A+t) ~ and the time development of x, 6 is governed by the second term 
in Eq. (ITOl) . At intermediate times tint with 1/|A+| ^ Unt ^ VI-^-l the velocities are 
approximately given by 













A--(C)ii 1 






\ (C)l2 / 



A_ - A_ 



Expanding R6/x in powers of {z — R)/R -C 1, we obtain R6/x = 1 + 0{{z — R)/R) 
which reflects the definition of rolling in mechanics. Thus, by the action of force and torque 
resulting from the tether bond, translational and angular velocity of the cell are adjusted and 
the cell starts rolling a short time after the first bond is formed. This can be seen in Fig. [5^, 
there the two velocities x, 6 are plotted as a function of time. Although Eq. ffTOj) correctly 
predicts the cell to arrest at large times, the linear approximation Eq. ([8]) is not applicable 
at these times, as can be seen from Fig. [51 In Fig. [5^ the cell velocities as predicted form the 
analytical solution to Eq. ^ are compared to the velocities which are obtained when Eq. ^ 
is numerically integrated. We see that for short times the linear approximation works quite 
well, until the cell has reached its rolling state (i.e., RQ/U ~ 1). At long time scales the 
approximate solution provides only a qualitative prediction of the stopping process. The 
figure shows that the higher order terms contributing to the bond force/torque lead to a 
much faster stopping than predicted by the approximation of these quantities with respect 
to first order in x, 9 (see Eq. (IHD). 

We now briefly discuss the case of a non-vanishing resting length Iq of the bond (a 
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FIG. 5: Comparison of the analytical solution for the stop dynamics mediated by a single bond 
Eq. (jlOp with values obtained from numerical integration of Eq. ^ {z = l.Oli?). All plotted 
quantities are in dimensionless form, i. e., force is expressed in terms of 6TTr]R'^^, torque in terms of 
QirrjE?^, velocity in 7i?, and time in I/7. For the dimensionless spring constant k = lO'^ was used. 
For a set of typical parameter values: R = 5 fim, 7 = 100 Hz, r/ = 10^^ Pas, k = 1 corresponds 
to 10^^ N/m. (a) Plot of the (angular) velocity x (6) as a function of time for three different 
calculations: from numerical integration of Eq. ([6]), from analytical derivative of Eq. (jlOp . and 
from numerical integration of Eq. ([6]) under the assumption that the initial bond length ro is the 
resting length of the bond lo = ro. For the time axis a logarithmic scale is used, (b) The same for 
the bond force (torque) Fb,x {TB,y)- 

non-vanishing resting length with /q ^ ^0 is used in our simulations as explained in the 
Appendix). For this we assume the height z — R oi the cell above the wall to be equal to 
the resting length of the tether bond. In that case the spring constant in Eq. ([7]) is replaced 
by K /t(l — /o/||rz — frll) (see Fig. H]). The leading order term in the power series of Fb^x 
is then of second order in x, 9. Thus the linear approximation of the equations of motion 
cannot be applied in this case and therefore we consider this case only numerically. The 
results are also shown in Fig. [5^. We note that the qualitative behavior is very similar to 
the case of zero bond resting length. Fig. [5^ shows that at short times the cell moves on 
with almost unchanged velocities, then the velocities adjust to R9 ~ x (for some time even 
R9 > ct) and finally the cell stops after about the same time as in the case with zero bond 
resting length. 

The linearized analysis of the rather complex motion of the cell under the action of a 
single bond might appear to be somehow simplifying, but it nicely reveals the mechanism 
that leads to cell rolling. This can best be understood from the time dependence of the 
torque that is exerted by the tether, see Fig. [Sjo. In our approximation the torque is given 
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FIG. 6: (a) Mean velocity {U), its standard deviation ajj, and mean angular velocity R{^}) as 
functions of the dimensionless rates 7r,eo. (b) (U), R{il,),au,aRQ as functions of the unstressed 
off-rate eo for two different on-rates vr. Parameters used for these simulations: R = 4.5 • 10~^ m, 
T = 293 K, 7 = 100 Hz, Ap = 0.05 • 10^ kg/m^, -q = 1.002 • lO^^ Pas, k = 1 • lO'^ N/m, 
ro = 1.0 • lO^^R , d = 5 ■ lO^^R, = 5000, Xc = 2 • 10^^^ m. The average was obtained over ten 
simulation runs of 20 s duration. 

by = K,{Rx — zR6) (see Eq. ([8])). Initially, x increases faster than R6 and a.s z ^ R the 
torque is positive, i. e., it supports the shear torque and the cell starts turning faster. At the 
same time the force RFb^x ~ ~TB,y slows down the translational motion. The maximum 
torque Tg = is reached when z9 = x, i. e., when the cell is approximately rolling. From 
that time on x and 6 increase at approximately the same speed and eventually the torque 
will become negative (if x < z9) and will act against the shear torque. Similar arguments 
hold when repeating the previous discussion with the exact expression for the torque given 
in Eq. ^. 

C. Cell motion at multiple bonds: classification of states of motion 

The previous analysis did only include a single bond which in addition was not allowed to 
rupture. In the presence of multiple bonds permanently forming and rupturing the situation 
is much more complex. Whether a cell is able to roll or not depends then on the one hand 
on external parameters like ligand-density, shear rate and viscosity. On the other hand it 
depends also on the internal parameters of the single receptor-ligand complex, which in our 
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state 


definition 


fypp m of ion 

rolling adhesion 
firm adhesion 
transient adhesion I 
transient adhesion II 


R{n)/{U) > 0.8 AND 0.95 > {U)/Um > 0.01 
{U) < 0.01 Um 

AND (Tu/(U) < 0.5 
0.01 < {U)/Um and R{Q)/{U) < 0.8 j^^^ ^^^j^J ^ 



TABLE II: Five stationary states of leukocyte motion. 

model are the on-rate kon, the off-rate ko and the detachment force F^. In the following we 
will present our results mainly as a function of the two internal rates in their dimensionless 
form, TT = kon/i and eo = k^/'^. In Fig. [6^ we plot the mean translational and angular 
velocities as well as the standard deviation of the translational velocity, au = a/ {U"^) — {U)"^ 
(where the average is an average over time and an ensemble of cells), in a large range of 
values for the dimensionless on- and off-rates. To further illustrate the dependence of the 
kinetic quantities on the on- and off-rate, in Fig. Eb the eo-dependence at fixed on-rate vr is 
re-plotted. At vr = 10, one nicely sees that with decreasing eg, translational- and angular 
velocities first approximate each other, i.e., U decreases and RVt increases. Then, both 
together decrease to zero at very low off-rates. At smaller on-rates (tt = 0.04), both U and 
RVL monotonically decrease with decreasing off-rate and RVt ^ U occurs only when both 
quantities are close to zero. The standard deviations of the velocities au-, (Trq are small 
for very low and very high off-rates. In between, they pass through a maximum which is 
located exactly at the transition from unperturbed motion to cell arrest {U ~ RQ ~ 0). We 
now summarize these qualitative observations by defining five different classes of stationary 
states of cell motion (see also Tab. HTll . 

Free motion: We call a cell to move freely if its speed is larger than 0.95 Uhd (an example 
is shown in Fig. 0(1)). Free motion as we define it does not imply that there are no bonds 
at all. The definition given by us rather allows also for bonds with a very fast dissociation 
rate (off-rate) or very small detachment forces. In this case existing bonds dissociate before 
they are stretched enough to apply forces that slow down the mean velocity of the cell below 
95 % of Um. 

Firm adhesion (arrest): This is the state when the mean translational velocity (U) is less 
than O.OlUhd (Fig. [Z(2)). This still allows for small jumps due to rare dissociation events. 
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Besides that tether bonds compensate shear force and torque (cf. Sec. IIII Bp . 

Rolling adhesion: The ratio R{Q) / (U) is larger than 0.8. As was shown in Sec. Ill Al this 
is well above the hydrodynamic maximum of this ratio in the limit z ^ R {i. e., when the 
cell touches the wall). Fig.[7](5) and Fig. [71^6) show two examples of computational leukocyte 
rolling. 

Transient adhesion: If none of these criteria applies we define the state as being transient. 
Within this category we distinguish two sub-classes according to the standard deviation au- 
By au/{U) < 0.5 the first sub-class {transient I) is defined, otherwise the cell's motion is 
in the sub-class transient II. 'Transient I' occurs if bonds form and rupture permanently, so 
that they reduce the (translational) velocity considerably below the hydrodynamic velocity. 
However, in this case the bonds do not last sufficiently long as to increase the ratio RQ/U 
above 0.8. Fig.[7](3) shows an example for this kind of motion. 'Transient II' is characterized 
by alternating periods of arrest and free motion which is illustrated in Fig. [3^4). 

As the kinetic quantities {U),{Q),(Tu vary continously with respect to n and eo, the 
classification given above is not unique. But it allows to clearly distinguish these states in 
an on- off' state diagram, i.e., the states are in general not degenerated. Other classifications 
of leukocyte states have been given before. For example, in the first paper on adhesive 
dynamics also five states of motion were defined ISj. In contrast to our definition, however, 
this classification was only qualitative. In computer simulations of adhesive dynamics or 
other numerical models for rolling adhesion, it is common to define states of motion like 
free, rolling and firm adhesion, but usually this is done like in an experimental context, 
which means that these classifications are based only on the mean translational velocity 
[31I . Is?! . A classification into immobile, rolling and detached and a corresponding state 
diagram has also been given by Bruinsma in a mean field __approach, which did not model in 
detail how force is distributed over the molecular bonds 



IV. STATE DIAGRAM OF LEUKOCYTE MOTION 

In order to determine the stationary state of motion for a given set of parameters, we 
repeatedly performed the following computer simulation. At a given set of parameter values 
we let the cell start at a height z = R + rQ and subsequently simulate its motion for 20 s. As 
a result of the downward acting buoyant force, which drives the cell even closer to the wall. 
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FIG. 7: Translational and angular velocities U and RQ, respectively, that give examples of the 
different states of leukocyte motion defined in this section. The labels (1-6) refer to the numbered 
points in the state diagram shown in Fig. [HI (1) Free motion, (2) firm adhesion, (3) 'transient I', 
(4) 'transient IF, (5,6) rolling adhesion. 



the wall ligands will be immediately within the capture range of the cell-receptors. Then, 



as shown in 



2l| . the mean time for receptor-ligand encounter is close to zero for typical 
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FIG. 8: (Color online) The on-ofF state diagram displays different stationary states of leukocyte 
motion obtained by simulations. The ordinate shows the dimensionless unstressed off-rate eo, the 
horizontal axis the dimensionless on-rate vr. Snapshots for the cell velocities at the (vr, eo) parameter 
values marked by numbered circles are shown in Fig. [71 Parameters used are the same as for Fig. [6l 

ligand and receptor densities found for leukocytes. Therefore, cell-wall interactions arising 
from bonds are assumed to influence the cell motion for the complete run of the simulation. 
To nevertheless rule out any initialization effects the mean values and variances for U, Q are 
only calculated for times greater than 4 s, i. e. 20% of the total run length of 20 s. To ensure 
proper classification of the state of motion, each simulation run is repeated at least ten times 
(each time with another randomly chosen receptor distribution) and each run contributes to 
the mean values (U), (Q) and their standard deviations (Ju,a^ (in the case that au is large 
even more than ten simulation runs were performed). The numerical time step was chosen 
to be between 10^^ — 5 ■ 10~^ (at a typical shear rate of 7 = 100 Hz, the lower limit of the 
numerical time step correspond to real time step of 10^^ s). The smaller range of time steps 
is chosen when high ligand densities or stiff bonds (large k) are considered (to avoid too 
large update steps at large direct forces in Eq. ([3])). Keeping the other parameters fixed and 
varying the rates vr, eo on a grid in double logarithmic scale, we can determine the different 
types of leukocyte motion in an on-off state diagram. 

Fig. [8] shows a first example of such a state diagram. The parameters used there are 
listed in the figure caption of Fig. [6] (as we keep the parameters R,Ta, Ap,ro fixed for all 
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the diagrams shown in this section, they are not exphcitly hsted for the following state 
diagrams) . 

All five states can be identified in Fig. [S]and in Fig. [7](l-6) example trajectories for each 
of these states are shown. In the limit of very large off-rates eg the cell moves freely, i. e., no 
matter how frequently bonds are formed, force cannot build up because dissociation occurs 
immediately. At very small off-rates the cell is in the state of firm adhesion, i. e., once a bond 
is formed it lasts long enough to stop the cell and to allow further bonds to be formed, thus 
stabilizing the cell in its rest position. In between these two limiting cases for the off-rate we 
find the other three states. From these the rolling state appears only for on-rates vr above 
a certain threshold. This confirms the conclusion drawn from early experiments with flow 
chambers that selectins are especially suited for rolling adhesion due to their fast on- and 
off-rates {cj. The state 'transient I' appears when the off-rate is too large to support rolling, 
but still too small to allow for free motion. When both on- and off-rates are small, the cell is 
in state 'transient 11'. In this state the cell stops most likely whenever a bond is formed due 
to the small off-rate. This results in periods of firm adhesion. The small on-rate, however, 
makes it rather unlikely that bonds are formed, which results in periods of free motion. 
Thus, the cell moves in a stick-slip like fashion. 

The dynamic state diagram Fig. [8] emphasizes that the molecular rates (the on-rate vr and 
the off-rate eg) are the main determinants of rolling adhesion. In addition, the state diagram 
also depends on the other parameters contained in the adhesive dynamics algorithm. In the 
following we discuss the qualitative dependence on some of these parameters. 

A. Medium viscosity 

The impact of viscosity is illustrated in Fig. [9K. There, the on-off diagrams for the 
fluid viscosities r] = 1, 3, 5 mPas are shown. The viscosity of the standard medium in which 
leukocytes are usually perfused through flow chambers is about 1 mPa s. For rj = 3 mPa s the 
states of motion are distinguished by areas filled with different grey scales. For rj = 1,5 mPas 
only the border lines between two states are shown. The figure legend explains which state 
can be found below a given line. For the sake of clarity, the original (12,13)-simulation grids 
for the three diagrams are omitted. Fig. [9ti clearly shows that the main effect of increasing 
the viscosity is the increase in the range where rolling is possible. More precisely, the larger 
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FIG. 9: (a) On-off state diagrams for different viscosities rj. The filled areas define the states of 
motion as determined from simulations at a reference viscosity r] = 3 mPas. The lines denote the 
borders between the states of motion at a lower and at a higher viscosity than the reference viscosity. 
The thin (thick) lines refer to a viscosity rj = 1 mPa s (5 mPa s) . The figure legend tells which state 
is found below the corresponding line, (b) States of motion at different spring constants k. The 
filled areas are the states of motion for the intermediate spring constant k = 5- 10~^ N/m. The thin 
(thick) lines refer to k = 5 • 10^^ N/m (k = 5 • 10^^ N/m). (Other parameters used for simulations 
in both (a) and (b): Number of receptors Ny. = 2500, ligand-ligand distance d = 2.5 • 10^^ R, in 
(a): K = 5 ■ 10^'^ N/m; in (b): rj = 1 mPas.) 



the viscosity, the lower the off-rate eo at which firm adhesion sets in. This effect results 
from the Bell model for bond dissociation. The shear stress 777 and the maximum force in a 
tether bond are proportional to the viscosity rj (the maximum force is the force which holds 
the cell at rest, see Fig. Hj). Therefore, an increase in viscosity from rj* to rj increases the 
off-rate like 

e = eo[exp(F^/F,)](^/^*), 



with the bond force at viscosity rj* . As a rough estimate for the bond force we use 
Fn ||Fs||/ cosx (see Fig. |4]for the definition of F5, x; for the angle we estimate x ~ 65° 
[gJ). Then, for the parameter values used here we have F^/Fd ^ 0.7 at r^* = 1 mPas. Thus, 



if at some viscosity rj firm adhesion occurs for off-rates smaller than a certain value e 



firm 
' 



I.e., 



for eo < eg**^™, then, we expect firm adhesion for rj* to exist at the same rate e. For rj* this 
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rate is estimated to be e^^^"^ exp(F^/ Fd)^^^* ^. With r]/r]* = 5, we therefore expect rolhng 
at T]* for eo > ef'^"'(r/*) ~ 156^*^(77 = 5t]*). A factor of roughly this order of magnitude 
between the off-rates at the border between rolhng and firm adhesion for 77 = 5 mPas and 
1]* = 1 mPa s can also be read off from Fig. [9^. On the other hand no essential shift in the 
borderline between the states rolling and the 'transient I' can be spotted when the viscosity 
is changed. This is the case as the referred line occurs at rather large eo values at which eg 
appears to dominate over the force (viscosity) dependent part in the Bell equation for the 
total off-rate e. These estimates suggest that the rolling state disappears at even smaller 
viscosities (one or two magnitudes smaller than that of water) than used in Fig. [9^. This 
is indeed observed for simulations in this viscosity range (data not shown). Note however 
that in this rage, the assumption of small Reynolds number might fail. 



B. Bond spring constant 

Fig. [9)d shows the state diagrams for three different spring constants i^suff, f^int, i^soft = 
5 ■ 10^^,5 • 10^'^, 5 ■ 10^^ N/m, respectively. The intermediate spring constant is of the 



same order of magnitude as the spring constant of the bond of P-selectin and its ligand [48] 



The softest spring constant Kgoft mimics the effect of soft microvilli |49l]. From Fig. [9b we first 
notice that the firm adhesion state for Hsuff occurs at smaller off-rates eo for small on-rates 
TT and at larger off-rates for large on-rates compared to the case of the intermediate spring 
constant hiint- A closer view identifies two effects that are responsible for this observation at 
small on-rates. First, the stiffness of the bond results in a small elongation which then leads 
to an obtuse bond- wall angle (the angle x Fig- HI)- The more obtuse this angle is, the more 
the bond must be loaded to compensate the shear force. In addition, the transport of the cell 
and thus also the bond extension is governed by the shear rate 7. A stiffer bond is therefore 



loaded faster. Both the faster loading 59|] and the larger bond force result in an effective 
increase of the off-rate e, which at small on-rates can be compensated by smaller off-rates 
eg. In contrast, at larger on-rates the fast dissociation is compensated by fast binding of 
new bonds and rebinding of just dissociated bonds. Therefore, the increase in the effective 
e plays a minor role. On the other hand the faster loading leads to a faster stop of the cell, 
so that to maintain rolling even larger rates eo are necessary in this range of vr. 

When we compare the intermediate spring constant with the soft one Kso/t, we see 
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FIG. 10: (a) On-off state diagrams for different numbers of receptors Nr. Nr = 1000 (thin lines) 
Nr = 2500 (filled areas) and Nr = 5000 (thick lines), ligand-ligand distance d = 5 • 10~^, k = 
10~^ N/m, 7] = 1 mPas. (b) State diagrams for different ligand-ligand distances d = 1 ■ 10^^ R 
(thin lines) d = 5 ■ 10~^ R (filled areas) and d = 2.5 • 10~^ R (thick lines), number of receptors 
Nr = 5000, K = 5 • 10^3 N/^^ ^ = X mPas. 



that the region of rolling shrinks for the soft case. In the case of the soft spring constant 
at larger vr rolling turns into 'transient I' already at smaller off-rates eo. In addition, rolling 
occurs still at much smaller on-rates than in the case of Kmt and i^suff-, respectively. Both 
of these observations have their origin in the larger tether elongation that is possible for soft 
bonds. This elongation effectively increases the contact area (i. e., the area on the cell surface 
that is less than the capture length tq away from the wall). More precisely, bonds can still 
only be formed within the contact region but can then exist also outside the contact region. 
If the off-rate eo is not too large, this effect leads to an effective increase in the number of 
available receptors, which explains the rolling in the region of smaller vr. For larger cq bonds 
rupture quickly also at small bond forces. In the case of soft bonds the bond force is then 
not sufficiently high as to reduce the translational velocity to induce rolling. The state of 
motion then rather appears to be 'transient I' or free motion. 
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C. Receptor number 



Before discussing the impact of the number of receptor patches on the cell surface, we 
first estimate the number of receptor patches within the contact zone given A^^- Supposing 
the cell touches the wall it immediately follows from basic geometrical considerations that 
on average Nr^contact = Nj. ro/2 receptor patches are in principle in capture range to wall 
ligands. On the other hand, the average number of ligands in range to receptors is calculated 
form the product of the projected area of the contact zone to the wall and the ligand density. 



The radius of the projected contact area is ~ y/2RrQ as tq <^ R. Thus the number of ligands 
in the contact zone is N^contact = 2nRrQ/cP. 

In Fig. [lOji we have Nr = 1000,2500,5000 and c/ = 5 ■ lO^^i?, i.e., in the contact zone 
Nr,contact ~ 5,12.5,25 and Ni^contact ~ 25 (for ro = 0.01/2). Therefore, the number of 
receptors limits the maximum number of bonds in all three cases. In addition, using as 
an estimate for the receptor patch density on the sphere pr = nr^N^/ [AttR^), we have 
Pr < 0.125, i.e., even for = 5000 receptor patches, they do not cover the contact zone 
completely. So, not every receptor necessarily encounters a ligand, and therefore the number 
of tether bonds is even less than Nr^contact- In fact, in the roUing state, we actually measure 
an average number of (1.3 - 2.6) for A^'^ = 1000, (2 - 5) for A^'^ = 2500, and (2 - 6) for 
Nr = 5000 existing bonds, respectively, depending on the actual on- and off-rate, which is 
less than the respective Nr^contact- Fig- flOb - shows that the larger A^^, the smaller on-rates 
TT are sufficient to support rolling. Not every receptor-ligand encounter (i.e., overlap of the 
receptor patch with a ligand) turns into a bond. This happens only with a probability 
depending on vr and the dwelling time of the encounter. Thus, the more encounter occur 
per time the more bonds will be formed at a given rate vr. As the encounter rate increases 
with increasing number of receptor patches 2l|, this also increases the average number of 
bonds. 

Fig. [TOk shows that the 'transient IF region expands with decreasing A^^- The main effect 
here is that for smaller A^^ the rolling turns into transient motion, while the border line 
between free and transient motion is much less effected by the decrease in the number of 
receptor patches. The large (vr, eo)-range for the transient states is a signature of few bonds 
being at work as we will see again when discussing the influence of ligand-ligand distance. 
Then, single tethers slow the cell down (depending on the off-rate they may either arrest 
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the cell some while, resulting in state 'transient IF, or just decelerate them resulting in 
state 'transient I'), but after dissociation it is unlikely that the current state of motion is 
supported by further bonds. However, as long as at least two bonds are possible this effect 
is partly compensated at large on-rates when the probability for receptor-ligand encounter 
to result in a bond is high. 

D. Ligand density 

To demonstrate the impact of ligand-ligand distance. Fig. [TUb shows the state diagrams 
for d = 10- 10^2 R, 5 ■ lO'^ i?, 2.5 • 10"^ R and A^^ = 5000. Using the expressions derived 
in the previous subsection we now have Ni^contact ~ 6, 25, 100, respectively, for the number 
of ligands in the contact zone. The mean number of receptors that may form a bonds is 
Ny^contact = 25. So, in principle Nr^contact limits the number of bonds in the two cases of higher 
ligand density and Ni^contact is the hmiting value only for the very low ligand density. When 
we measure the average number of bonds at high on-rates and relatively small off-rates (i. e. 
in the rolling and firm adhesion region) we find for the intermediate ligand density with 
d = 5 - 10~^ R about five bonds, and for the large ligand density with d = 2.5 ■ 10~^ R about 
fourteen bonds on average. The later value is slightly larger than N^^contact = 12.5. This can 
be explained with the elasticity of the bonds that — once formed — allows them to also exist 
beyond the contact zone. Thus only at very high ligand densities with Ni^contact ^ Nr^contact, 
the number of receptors limits the maximum number of bonds. For ligand densities with 
Ni^contact ~ Nr^contact & dccrcasc in d still leads to an increase in the average number of bonds 
as this increases the rate of receptor-ligand encounters. 

The basic effect of decreasing the ligand density as illustrated by Fig. [TOb is the shift of 
all state from the upper left towards the lower right in the on-off state diagram. For example 
rolling is only supported at larger on-rates and smaller off-rates when the ligand density is 
decreased (i.e., d is increased). Also the border line between the state of free and transient 
motion is notedly shifted between the two extreme cases of ligand density simulated for 
Fig. [TOb . In fact, the shift of this border line is much more pronounced than in the previous 
discussed case where the number of receptor patches was reduced. The simple reason for 
that is that in Fig. [TUk Nr^contact is changed by a factor of five whereas in Fig. [TUb the 
Ni^contact is changed by a factor of almost twenty. 
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FIG. 11: Determination of the dynamic on-rate vr from the determination of the state of motion. 
At the fixed off-rate eg = 0.16, the ratio R{Q)/{U) is plotted as a function of the on-rate vr for 
different ligand-to-ligand distances d ranging from d = I - 10-2 to d = 5 • 10^. In the vr- range 
above the solid line at R{il.)/{U) = 0.8 the ceU is rolling. In the inset the on-rate vr* at which 
R{Q) I (U) = 0.8 is plotted as a function of d. (Other parameters used: = 2500, d = 2.5 • 10"^ R, 
7 = 100 Hz, Xc = 2 • 10"" m; K = 5 • 10"^ N/m; r/ = 1 mPas.) 



Similarly to the case of reduced number of receptor patches, a small ligand density results 
in an increased (tt, eo) range for transient adhesion. Thus, we note again that at a reduced 
rate of receptor-ligand encounters, rolling tends to be converted to transient motion. We 
also note that if several parameters are changed at the same time, the overall effect can 
be qualitatively understood by superimposing the effects of the single parameter changes 
described above (data not shown). 



E. Application: experimental determination of the on-rate 

The on-off state diagrams introduced here display the state of motion of a cell or a 
receptor covered micro-bead for a given set of experimental parameters in a wide range of 
possible on- and off-rates of the receptor-ligand pair. Because the molecular properties of 
a given receptor-ligand pair can be assumed to be given a priori, each pair corresponds to 
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one point in the on-off state diagram. In the following we show how the considerations 
regarding the rolling state presented here can be used to determine the dynamic on-rate 
vr of a receptor-ligand pair by changing the external parameter. For this line of reasoning 
the off-rate eo is assumed to be known, e.g. from dynamic force spectroscopy experiments. 
Fig. [TT] shows the ratio of the mean velocities R{Q) / (U) of the cell as a function of vr for 
different ligand-to-ligand distances d. We see that for a given d the cell is slipping if the 
on-rate is small. Increasing the on-rate turns the cell's motion into rolling, i.e., at some 
critical on-rate vr* the ratio R{Q)/{U) reaches the value 0.8 which we defined for rolling. 
The smaller (larger) the ligand density (ligand-to-ligand distance d) the larger is the critical 
on-rate vr*. More quantitatively this is shown in the inset of Fig. [TT] where vr* is plotted as a 
function of d. If one now changes the ligand density in a series of flow chamber experiments 
from low to high values and by doing so determines the ci-value at which the motion of the 
cell turns from transient to rolling motion, then one can read of from the inset of Fig. [11] 
the on-rate vr of the receptor-ligand bond. 

The described procedure does not work for very small on-rates for which rolling is not 
possible even for very large ligand densities (i. e., d tq). On the other hand one can always 
find a ligand distance at which rolling is not possible whatever value the on-rate takes. At 
very low ligand densities d ^ ro the reaction is limited by the encounter rate and even for 
very large rates vr rolling is not possible. 

V. DISCUSSION 

In this paper we introduced a new version of the adhesive dynamics algorithm. In con- 
trast to earlier work, our approach also includes the diffusive motion of cells resulting from 
thermal fluctuations. This allowed us to spatially resolve receptors and ligands. An imme- 
diate advantage of this approach is that the single bond on-rate kon can be chosen to be 
independent of the relative motion of cell and substrate. In this work, we focused on the 
different dynamic states of motion which can be identified on the cellular level. We first 
noticed that the action of a single bond not only slows the cell down, but also changes the 
motion from slipping (which is the case for cells moving free in hydrodynamic flow) to rolling 
in the sense RVt/U ^ 1. In the case of multiple bonds, rolling can also be observed regarding 
the mean values of the velocities, i. e., R{VL) / (U) ^ 1 at proper rates of bond formation and 
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rupturing. This motivated us to define the state of rolhng adhesion as R{fl) / (U) > 0.8. By- 
extending these observations we defined five distinct states of stationary cell motion. These 
states were then displayed in so-called on-off state diagrams which showed the impact of 
different molecular rates. In addition, we investigated the effect of external parameters. For 
example, we showed that the cellular motion is changed considerably when the viscosity of 
the medium is changed. 

Our work shows that different dynamic states of cell motion can be defined in a systematic 
and quantitative way. In particular, calculations of state diagrams allow us to obtain a 
complete understanding of the way molecular and other parameters determine motion on 
the scale of a cell. Similar approaches have been taken before 57|, |58| , but without a proper 
definition of rolling in the mechanical sense. Our computer work now shows that including 
the rotational degrees of freedom allows us to investigate the transition from hydrodynamic 
slipping to bond-mediated rolling in a more detailed way. Experimentally it is certainly 
a challenge to obtain similar data in fiow chamber experiments for micron-sized beads or 
cells. In principle, making use of recent nanotechnological developments one could attach 
receptors to micron-sized beads which are covered with anisotropic surface layers 6C 



m. If 



these layers are anisotropically refiective, rotational motion of the spheres can be recorded. 
For cells, one would have to track surface or intracellular markers (e.g. mitochondria or 
nuclei). 

We also suggested a new procedure to experimentally determine values for on-rates from 
monitoring cell motion (cf. Sec. IIVEI) . In contrast to the off-rate, which can be measured 
e.g. from dynamic force spectroscopy experiments, it is very difficult to infer values for the 
on-rate kon in cell adhesion, where both partners have to be attached to appropriate surfaces. 
Here, also one could make use of recent developments in nanotechnology. As we have shown 
above, the ligand-to-ligand distance d is a crucial control parameter in our system. Recently, 



is has become possib 
arrays of gold dots 



e to control this parameter using nanopatterend and biofunctionalized 
62I . Therefore, in the future our predictions regarding the effect of 
ligand positioning might be compared directly to experimental data, especially if compared 
with a measurement of the rotational degrees of freedom. 

For conceptual and computational simplicity, here we have used the Bell model for the 
force dependence of the off-rate ^o//- In principle, it would be easy to include more com- 
plicated rupture scenarios, like the catch-slip behavior recently reported for both P- and 
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L-selectin 



BE3. 



There is good reason to believe that this molecular behavior is essential 



for the physiological function of these molecu 
suggested to explain this behavior 
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es and different theoretical models have been 



65| . If combined with our modeling framework 



for adhesive dynamics, these models might be tested against experimental data from flow 
chamber experiments. Such an approach would have the big advantage that it avoids testing 
single molecules outside their physiological environment, which is especially problematic for 
adhesion receptors which usually are embedded in the plasma membrane and regulated by 
the cytoskeleton. 

Further possible extensions of our simulation framework include models for cell deforma- 
bility and hydrodynamic interactions between cells. Cell deformations in free flow should 
become relevant only at shear rates well above 100 Hz 22|. In adhesion, cell deformation 
depends also on the number and strength of adhesion bonds. In the case of strong adhe- 
sion, it is well known that also viscoelastic changes occur, including elongation of microvilli 
491]. For computational and conceptual simplicity, here we have focused on the case of 
moderate shear flow and adhesion, when the effect of cell deformability is small. Moreover 
deformations are irrelevant for rigid microbeads, which have been shown to results in similar 



physical effects as described for cellular systems 111, Q, ll 



In order to combine elastic 



effects and hydrodynamics, a very powerful framework is provided by multi-particle collision 



dynamics 



661], which in principle also would allow to include non-laminar flow conditions 



and hydrodynamic interactions between cells. However, these effects can be safely avoided 
in flow chamber experiments by using sufficiently small shear rates and cell numbers. As 
explained above, an elegant way to test our theoretical predictions experimentally would be 
the combination of appropriately coated microbeads with nano-structured substrates. 

Finally we comment on the applicability of our approach to other systems that are based 
on the stochastic interplay between transport and adhesion. In the work presented here, we 
have considered a situation where receptor and ligands are both tethered to macroscopic 
surfaces. However, similar physical processes are relevant if the two molecular binding 
partners are free in solution, for example in affinity measurements 6?| or shear-induced 
adhesion of blood-clotting factors [68]. Another biological system for which the concepts 
for bond formation, bond rupture and transport discussed here can be applied is the cargo 



transport by multiple molecular motors 
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701 ]. Compared with rolling adhesion, the cell 



is replaced by the cargo (e.g., a vesicle), the cell-anchored receptors by molecular motors 
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and the substrate- anchored hgands by binding sites on the filament. In contrast to the case 
of rolhng adhesion, now the system is not driven by some external force, but the molecular 
motors actively step forward and pull their cargo against an external viscous friction force. 
Our approach can also be applied to non-biological systems. For example, it is well known 
that erratic motion also occurs in the context of sliding friction. A flat slider that is pulled 
above a plain wall exhibits stick-slip motion in some range of the pulling velocity. It was 
shown recently that this special kind of motion can be explained assuming the stochastic 
formation and rupture of molecular bonds between the slider and the wall TlJ. Thus the 
stick-slip motion of sliding friction and the erratic movement of cells in rolling adhesions 
seem to be based on the similar physical principles. 
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APPENDIX A: BOND DYNAMICS ALGORITHM 

The sphere's motion is described by Eq. ([3]). If no bond between the sphere receptors 
and wall ligands exists, we take for the direct force F-^ only gravity into account, i.e., the 
six-dimensional force and torque vector is given by F^ = {—AmgezjO) with g the earth 
acceleration constant. Receptor-ligand bonds lead to additional contributions to both the 
force- and momentum-part of F^. More precisely, a bond between a ligand located at r/ 
and a receptor located at r,, on the sphere's surface (see Fig. [1]) pulls with a force 

FB = hF{n), h-=7r r,, := Hr^ - r^||. (Al) 

III"/ — Trll 

F{x) is the force extension curve that describes by what force the bond must be pulled to 
stretch it up to a total length x. Here, we consider the bonds to be semi-harmonic springs 
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(cable model) 

= k{x - lo)e{x - /o), e(x) :=<("' ^ V (A2) 

else 

with /o the resting length and k the spring constant. The cable model is the simplest model 
for polymeric tethers. In the cable model a bond behaves as a spring only if it is stretched 
(extension larger than the resting length), otherwise the bond exerts no force on the sphere. 



Treating t 
regime 



le receptor-ligand complex as a harmonic spring works fine in the small extension 
48j. For large extensions the force extension curve for polymers is supposed to grow 
much faster than linear, and when the bond extension approaches the total contour length 
of the receptor-ligand complex it even diverges (strain stiffening). However, typical bonds 
are weak and their rupture probability increases exponentially with force. Therefore, we 
expect bond extensions to be restricted to the linear regime. As the bond force pulls on the 
sphere's surface also a torque results 

Tb = r X FB(rfe), 

where r is the connection vector from the center of the sphere to the point on its surface 
where the receptor is attached (see Fig. S]). Thus, the total force and torque contribution to 
F"^ by the bonds is 



i=l 



^g.F(r^)(rLFxra, (A3) 



with Nr the total number of receptors and = 1 if the ith receptor forms a bond and zero 
otherwise. The qi,i = l,...,Nr are stochastic variables. Thus, the contribution Eq. (1A3P 
lets the direct force F^ also become a stochastic variable. 

With this at hand we can now define the adhesive dynamics rules, applied in each update 
step At (cf. Fig.[T]): 

(i) The sphere's position and orientation is updated according to Eq. ([3]) (for an explicit 



description see Ref. 22l|). 



(ii) The receptor positions in the flow chamber coordinate system are calculated. 
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(iii) Each inactive receptor is represented by a capture ball with radius ro ^ 1. 

(iv) If the distance between a receptor and any ligand is < tq a bond is established with 
probability = 1 — exp(— At • kon), then the resting length of the bond is set to the 
receptor-ligand distance at the instance of bond-formation (i.e., the bond force at the 
moment of bond formation is zero) and is stored together with the ligand position. A 
bond can only be formed if the corresponding receptor and ligand are not already part 
of another bond. 

(v) For each active bond, the contribution to F-^ is calculated. 

(vi) Each existing bond dissociates with a rate given by the Bell equation, Eq. ([5]). Thus, 
each bond ruptures with probability poff = 1 — exp(— At ■ koff{F)), where F is the 
instantaneous force acting along this bond. 

When a bond has ruptured, both the receptor and the ligand can form a new bond in the next 
time step according to the rule (iv). As for the resting length Iq of a bond always Zq < ''"o ^ 1 
is true, modeling bonds as harmonic springs in both the extension and compression regime 
would not make much difference to the results that are obtained by the cable model. Given 
the probability for bond formation or rupture po„ or Poff, respectively, a standard Monte- 
Carlo technique is used to decide whether the action happens or not: Using a pseudo-random 
number generator a random number rand from the uniform distribution in the interval [0, 1] 
is drawn. If then Pon/off > rand the respective action takes place, otherwise not. 
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